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Abstract. In the framework of the renormalization-group theory of critical 
phenomena, a quantitative description of many continuous phase transitions can be 
obtained by considering an effective theories, having an A^-component fundamental 
field and containing up to fourth-order powers of the field components. Their 
renormalization-group flow is usually characterized by several fixed points. We give 
here strong arguments in favour of the following conjecture: the stable fixed point 
corresponds to the fastest decay of correlations, that is, is the one with the largest 
values of the critical exponent 77 describing the power-law decay of the two-point 
function at criticality. We prove this conjecture in the framework of the e-expansion. 
Then, we discuss its validity beyond the e-expansion. We present several lower- 
dimensional cases, mostly three-dimensional, which support the conjecture. We have 
been unable to find a counterexample. 
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1. Introduction 

In the framework of the renormahzat ion-group (RG) approach to critical phenomena, 
a quantitative description of many continuous phase transitions can be obtained by 
considering effective Landau-Ginzburg-Wilson (LOW) theories, having an A^- 
component fundamental field $j and containing up to fourth-order powers of the field 
components. The fourth-degree polynomial form of the potential depends on the 
symmetry of the system. LGW theories generally present several fixed points (FP's) 
which are connected by RG trajectories. See, for example. Fig. [1], which shows the 
RG flow in the example of four fixed points. Among them, the infrared stable FP 
determines the asymptotic critical behaviour of the corresponding statistical systems. 
FPs are determined by the zeros of /3 functions. The stability of the FP is related to the 
eigenvalues of its stability matrix: if all eigenvalues have a positive real part, then the 
FP is stable. An interesting question is whether a physical quantity exists such that the 
comparison of its values at the FPs identifies the most stable FP. In two dimensions, the 
central charge is such a quantity: the stable FP in unitary theories is the one with the 
least value of the central charge [1]. But, despite several attempts and some progress, 
see, for example, Refs. [2], no conclusive results on the extension of this theorem to 
higher dimensions has been obtained yet. 

In this paper we give strong arguments in favour of the following conjecture: 
In general theories with a single quadratic invariant, the infrared stable FP is 
the one that corresponds to the fastest decay of correlations. 

Therefore, it is the FP with the largest value of the critical exponent r] which 
characterizes the power-law decay of the two-point correlation function W^'^\x) at 
criticality, 

W'^'i-)^^- (1) 

The exponent rj is related to the RG dimension of the field, (i$ = (rf — 2 + ri)/2. 

The conjecture holds in the case of the 0(A^)-symmetric theory. Indeed, below 
four dimensions, the Gaussian FP, for which 77 = 0, is unstable against the non-trivial 
Wilson-Fisher FP for which 77 > 0. We recall that the positivity of t] in unitary theories 
follows rigorously from the spectral representation of the two-point function [5]. 

In the absence of a sufficiently large symmetry restricting the form of the 
potential, many quartic couplings must be introduced — see, for example, Refs. [U HJ 
El El El E]- The Hamiltonian of a general $^ theory for a A^-component field can be 
written as 

i i ' ijkl 

The number of independent parameters and Uijki depends on the symmetry group 
of the theory. An interesting class of models are those in which J2i is the unique 
quadratic polynomial invariant under the symmetry group of the theory. In this case, 
all are equal, = r, and Uijki must be such not to generate other quadratic invariant 
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terms under RG transformations, for example, it must satisfy the trace condition [7] 
J2iUiiki oc 6ki. In these models, criticality is driven by tuning the single parameter 
r, which physically may correspond to the temperature. All field components become 
critical simultaneously and the two-point function in the disordered phase is diagonal, 
that is. 



These $^ theories have generally several FPs. Our conjecture applies to such class of 
unitary models. We do not consider nonunitary limits, such as ^ 0, which are 
relevant to describe the critical properties of spin systems in the presence of quenched 
disorder. 

Actually, although within the e-expansion around d = 4 one can prove that only one 
stable FP exists, in ci < 4 general theories may have more than one stable FP with 
separate attraction domains. The rj conjecture should be then refined by comparing FP 
that are connected by RG trajectories starting from the Gaussian FP: among them the 
stable FP is the one with the largest value of t]. 

It was already observed in Ref. [Ij that, within the e-expansion, the stable FP of 
theories with four FPs, is the one with the largest value of 7]. Here, we extend this 
e-expansion result to an arbitrary $^ theory. Then, we discuss the validity of the r) 
conjecture at fixed dimension d < A, where it remains a conjecture. We present several 
lower-dimensional cases, mostly three-dimensional, which support the conjecture. We 
have been unable to find an analytical or numerical counterexample. 

Finally, we extend the conjecture to multicritical points in models with several 
independent correlation lengths (different r^) that diverge simultaneously. In this 
situation, the exponent i] is replaced by a matrix and the conjecture applies to the 
trace of the matrix. However, the empirical evidence, beyond the e-expansion is more 
limited. 

The paper is organized as follows. In Sec. [2|, we prove the rj conjecture within 
the e-expansion for the most general <l>^ with a single quadratic invariant. In Sec. [Sj 
we discuss several lower dimensional examples, showing that in all cases the stable FP 
is the one with the largest value of rj. In Sec. HI we discuss the extension of the rj 
conjecture to $^ theories describing multicritical behaviours. In Sec. [5] we discuss the 
Gross-Neveu- Yukawa model, and show that also in this case the infrared stable FP is 
the one with the fastest decay of the critical two-point function of the boson field, within 
the e-expansion and in the large- A^/ limit for any dimension. 

2. Proof within the e-expansion 

Expansions in powers of £ = 4 — can be most easily obtained within the minimal- 
subtraction scheme 0, where the RG functions are computed from the divergent part 
of correlation functions [3], 



W^fix -y) = ($,(x)<l>,(i/)) = 5,,W^^\x - y). 



(3) 



Pijkl{gabcd) = 



Vidabcd) = fJ' 



dlnZ^ 



djji 



(4) 
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where gijM are the renormahzed couphngs corresponding to the quartic parameters Uijki- 
We consider only Hamiltonians that have a symmetry such that the quadratic 
invariant in the field is unique and the two-point function in the disordered phase thus 
diagonal. As a consequence of the diagonal property of the two-point correlation function 
in the disordered phase, the tensor Uijki has special properties that take the form of 
successive constraints in the perturbative expansion. At leading order one finds [7] 

1 

Pijkl(,9abcd) ^9ijkl ~l~ ^ ^ ^ (^9ijmn9mnkl 9ikmn9mnjl 9ilmn9mnkj) (5) 

±U/I rn,n 

The RG function associated with the field dimension can be inferred from the function 

V{9ijkl) = gjY(^4^-)4 X! 9ijkl9ijkl , (6) 



this form resulting from the diagonality condition 

^ ^ 9iklm9jklm ^ij ^ ^ 9klmn9klmn ■ (7) 
k,l,m k,l,m,n 

One can easily verify that the general expression of the one-loop /^-function derives 
from a potential [10]. Indeed, 

n ( X _ dU{9abcd) /„v 

Pijkl[9abcd) — ^ , [O) 

09ijkl 

U{9abcd) = -| E 9^,kl9^,kl + E 9^,kl9kln.n9n.n^, 

i,j,k,l ^ ' i,j,k,l,m,n 

Such a property, which has been verified to two-loop order, is shared, at leading order, 
by other field theories and would deserve a more systematic investigation. 

A detailed discussion of the properties of the RG fiow within the e-expansion can 
be found in Ref. [11]. Here, we list a number of consequences of Eq. ([8]) within the 
e-expansion. Note that none of these properties depends on the condition ([3]). 

(i) The potential decreases along a RG trajectory and thus FPs are extrema of the 
potential. In particular, if two FPs are (asymptotically) connected by a RG trajectory, 
the stablest FP corresponds to the lowest value of the potential. 

(ii) The eigenvalues of the matrix of first order partial derivatives of the (3 functions 
(stability matrix) at a FP are real. 

(iii) Stable fixed points are local minima of the potential, that is, the matrix of 
second derivatives of U{9) is positive. 

Moreover, two additional properties depend of the special cubic form of the one-loop 
potential ([H]) (we give the proof in the appendix): 

(iv) There exists at most one stable FP. 

(v) The stable FP corresponds to the lowest value of the potential U{9). 

The latter properties are not necessarily valid beyond the e-expansion. For example, 
in the physical dimensions d = 3,2 LGW $^ theories may have more than one stable 
FPs with separate attraction domains. This possibility does not occur within the e- 
expansion, that is, close to 4D, but it does not contradict general RG arguments and it 
is found in some cases, in particular when different regions of the Hamiltonian quartic 
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parameters are related to different symmetry breaking patterns. Next section we sfiall 
mention a few examples where this occurs. 

In the framework of the ^-expansion, we now show that the stable FP (or at least 
the most stable one) corresponds to the largest value of the exponent ri, and thus to the 
case where correlation function has the fastest decay at large distance. For any fixed 
point g^jki, the equations 

PiM9*a,J - = (9) 

implies 

3 

9ijkl9klmn9mnij 

(10) 

i,j,k,l V / i,j,k,l,m,n 

and, thus, at leading order, 

U{9*abcd) = E 9ijkl9ijkl (11) 

i^jyk^l 

which is negative, thus lower than the Gaussian FP value. At leading order, the exponent 
77 is then given by 

11 11 
^ " 6iV (47r)4 ^^3i3ki9ljki = ~iV£ [AtxY^^^*"^"'^^' ^^^^ 

As we have shown, the stable fixed point corresponds to the lowest value of U . It 
thus corresponds also to the largest value of the exponent rj: therefore, the correlation 
functions corresponding to the stable fixed point have the fastest large distance decay. 

The validity of this result beyond the e-expansion remains a conjecture. In next 
section, we discuss several checks in lower dimensions, mostly in three dimensions. 

3. Several verifications of the r] conjecture 

In this section, we discuss the RG flow, in c? < 4 dimensions, of several <l>*^ theories with 
a single quadratic invariant, but more than one quartic terms. As we shall see, in all 
examples considered below, the stable FP of the RG flow, within regions connected by 
RG trajectories starting from the Gaussian FP, is the one corresponding to the largest 
value of 77. 

3.1. The 0{M) ® 0{N) model 

We first consider the 0{M) ®0{N) <l>^ model corresponding to the hamiltonian density 

T^m^ai)' + r$L] + f E ^l] (13) 

a,i ' \o,,i J 

a,i,b,j 

where is a M x N real matrix (a = 1, M and i = 1, N). The symmetry of this 
model is 0{M)(^0{N). 
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Figure 1. RG flow of the 0{M) (g) 0{N) model in the large-iV limit. 



3.1.1. The large-N limit The 0{M) 0{N) $^ model can be solved in the large- iV 
limit for any fixed M p!2l [T3] . For any M > 2, one finds four FP's: the Gaussian FP, 
the Heisenberg 0{M x A^) FP, and two new FP's which we call chiral (C) and antichiral 
(A). Fig. [1] shows a sketch of the RG fiow in the quartic-coupling space. In the large- A^ 
limit, for any M > 2 and for any 2 < d < 4, the stable FP is the chiral one, and all 
other FP's are unstable. The large-A^ critical exponent t] has been calculated for all 
FP's: 



where 



and 



V 



ed 



N 



+ 0{l/N' 



AT{d 



r(2 - rf/2)r(rf/2 - i)r(rf/2 - 2)r(rf/2 + 1) 



(14) 



(15) 



(16) 



Gaussian 

1/M Heisenberg 

(M + l)/2 chiral 

(M - 1) (M + 2) /(2M) antichiral 

0(1/A^^) calculations can be found in Ref. [13]. One easily verifies that for any M and 
2 < d < A the value of r] at the stable chiral FP is the largest. 



3.1.2. d = 3 results from high-order FT perturbative analyses Several results have 
also been obtained concerning the RG fiow of the three-dimensional 0{M) ® 0{N) 
models at finite values of M, A^. Below we restrict ourselves to the case M = 2, that is, 
to 0{2)<^0{N) models. The cases A^ = 2, 3 are physically interesting because they could 
describe transitions in non-coUinear frustrated magnets, the superfiuid transition in ^He, 
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etc.... See, for example, Refs. [THE] and references therein. We have to distinguish the 
cases f > and vq < 0, because they lead to different symmetry breaking patterns: 



The M-axis plays the role of a separatrix and thus the RG flow corresponding to f o > 
cannot cross the u-axis. The relevant FPs of models with the symmetry-breaking pattern 
(JTTj) lie in the region v > 0, where v is the renormalized quartic coupling associated with 
Vq, while the relevant FPs of models with the symmetry-breaking pattern ffTSj) lie in the 
region f < 0. 

All 0(2) (g) 0{N) models contain the Gaussian FP and the Heisenberg 0{2N) FP. 
They are both unstable. The relevant perturbation at the 0{2N) FP, which makes it 
unstable for any > 2, is related to the w-term in the Hamiltonian, which is a particular 
combination of quartic operators transforming as the spin-0 and spin-4 representations 
of the 0(2A^) group. Any spin-4 quartic perturbation is relevant at the 0{K) FP for 
K > 3, since its RG dimension 2/44 is positive for A' > 3 [15]. Therefore the 0(2A^) 
FP is always unstable in the RG flow of the 0(2) (S) 0{N) models (actually this result 
extends to any 0(M) ® 0{N) model with M > 2). In particular, 7/4^4 ^ 0.11 at the 
0(4) FP and 7/4,4 ~ 0.27 at the 0(6) FP. 

The RG flow of three-dimensional 0(2) (g) 0{N) models has been investigated 
by computing and analyzing high-order perturbative series within the massive zero- 
momentum (MZM) and massless MS schemes, respectively to six and flve loops 
[T6| [T71 [T8| [T9] . Some results are reported below. 

RG flow for N = 2 and v < 0. One flnds a stable FP in the region f < 0, which 
is in the XY universality class [HI E] • We recall that the other relevant FPs are the 
Gaussian FP and the 0(4) FP. The best available estimates of rj for 0{N) models are 
reported in Tabled! They support the r] conjecture, which would require rjxy > '7o(4)- 
Indeed the best available estimate are rjxy = 0.0381(2) and ?7o(4) = 0.0365(10). 

RG flow for N = 2 and v > 0. The analysis of the high-order MZM and MS 
expansions provide a rather robust evidence of the existence of another stable chiral FP 
for V > [T6l [T7| [19]. This has been conflrmed by MC simulations of a lattice <I>^ model 
[in]. This FP is not connected with the one found in the f < region, because the line 
V = is a separatrix. The estimates of rj at this stable FP are: rjch = 0.09(1) from MZM 
and rjch = 0.09(4) from MS. These results must be compared with the values of rj at the 
other FPs connected by RG trajectories in the region f > 0, which are the Gaussian 
and the 0(4) FPs. Again the conjecture is verifled because rj^h > Vo{4,) > 0. 

RG flow for N = 3 and f < 0. There is a stable FP with attraction domain in the 
region f < [I8]. The corresponding estimates of r] are: rj = 0.079(7) from MZM and 
?7 = 0.086(24) from MS. The other FPs are the Gaussian FP and 0(6) FP, which have 
much smaller values of r], in particular ?7o(6) = 0.029(3) [28] . 

RG flow for N = 3 and v > 0. There is a stable FP also in the region f > |19j . 
The corresponding estimates of t] are: 77 = 0.10(1) from MZM and r] = 0.08(3) from 



0(2) (g)O(A^) ^ 0(2) (g)0(A^-2) for i;o > , 
0{2)®0{N) ^ Z2®0{N -1) fort;o<0. 



(17) 
(18) 
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Table 1. Best theoretical estimates of the critical exponents for three-dimensional 
0{N) models. Concerning the methods, IHT indicates high-temperature expansion 
of improved lattice models with suppressed leading scaling corrections; MC indicates 
Monte Carlo simulations; FT indicates field-theoretical methods based on perturbative 
expansions; MC-I-IHT exploits a synergy of IHT and MC. Other results can be found 
in Refs. [31[B] 



N 


V 




Method 


Ref. 




1 


0.63012(16) 


0.0364(2) 


IHT 


[20] 






0.63020(12) 


0.0368(2) 


MC 


EH 






0.6304(13) 


0.034(3) 


FT 


m 


[23] 


2 


0.6717(1) 


0.0381(2) 


MC+IHT 


[24] 






0.6703(15) 


0.035(3) 


FT 


[21 


[23] 


3 


0.7112(5) 


0.0375(5) 


MC+IHT 


[25] 






0.7073(35) 


0.0355(25) 


FT 


m 


[23] 


4 


0.749(2) 


0.0365(10) 


MC 


[26] 






0.741(6) 


0.0350(45) 


FT 


[22] 




5 


0.779(3) 


0.034(1) 


MC 


m 






0.764(4) 


0.031(3) 


FT 


[28] 




6 


0.789(5) 


0.029(3) 


FT 


[28] 




8 


0.830 


0.027 


FT 


[29] 




large N 


1 - 32/(37r2iV) 


8/(37r2Ar) 


1 /TV exp 


m 





MS. This values are again much larger than the values of r] of the unstable Gaussian 
and 0(6) FPs. 

Note that the stable three-dimensional FPs of 0(2) ® 0{N) models with = 2, 3 
do not exist close to four dimensions, apart from the one for N = 2 and v < 0, see, for 
example, Ref. [H]. Thus, they provide a rather non-trivial check of the rj conjecture, 
because the three-dimensional RG flow differs qualitatively from the RG flow close to 
d = 4, which is obtained from the e-expansion. 

3.2. $^ theory with cubic anisotropy 

This theory is relevant for magnets with cubic anisotropy. Its Hamiltonian density is 

^ = ^ E [{dM' + r<l>l] + (e ' + 1, -0 E 

where <l>j is a A^-component field. The RG flow of this Hamiltonian has four FP's: the 
trivial Gaussian one, the Ising one in which the N components of the field decouple, 
the 0(A^)-symmetric and the cubic FP's. The Gaussian FP is always unstable, and so 
is the Ising FP for any number of components A^. Indeed, at the Ising FP one may 
interpret the cubic Hamiltonian as the Hamiltonian of Ising systems coupled by the 
0(A^)-symmetric interaction. The coupling term / d'^x^^^^ with i j scales as the 
integral of the product of two operators ^f. Since the <I>^ operator has RG dimension 
1/z// — indeed, it is associated with the temperature — the combined operator has RG 
dimension 2/uj — d = aj/uj, and therefore the associated crossover exponent is given 
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Figure 2. RG flow of the cubic model. A'^c < 3 in d = 3, see, for example, Ref. [31] . 

by = a/, independently of A^. Since a/ > 0, the Ising FP is unstable independently of 
A^. On the other hand, the stability properties of the 0(iV)-symmetric and of the cubic 
FP's depend on A^. For sufficiently small values of N, N < Nc, the 0(A^)-symmetric 
FP is stable and the cubic one is unstable. For > Nc, the opposite is true: the RG 
flow is driven toward the cubic FP, which now describes the generic critical behaviour 
of the system. Fig. [2] sketches the flow diagram in the two cases N < Nc and A^ > Nc- 
In d = 3 Nc ^ 2.9, see for example, [311 [6] and references therein. Therefore, the 0{N) 
FP is stable only for N = 2. Let us now examine the various cases in more detail. 

The case N = 2. We have four FPs: the Gaussian FP, the Ising FP along the 
u = axis, the cubic FP for t> < which turns out to be equivalent to an Ising FP, and 
the 0(2) FP which is the stable one. Again the rj conjecture does not fail: it requires 
VxY > "/Rising and this is verifled by the best three-dimensional estimates reported in 
Tabled! i.e. rjxY = 0.0381(2) and T^ismg = 0.0364(2). 

It is also worth mentioning the two-dimensional RG flow of this theory, where there 
is a line of stable flxed points connecting the Ising and XY FPs [32l [33] , with central 
charge c = 1 (in this case the central charge associated with the Ising FP is c = 1 
because it represents two decoupled Ising models). Along this line the critical exponent 
7] does not vary, t] = 1/4, but the correlation-length critical exponent goes from u = 1 
to z/ = oo. 

The case N = 3. We have four FPs: the Gaussian FP, the Ising FP along the u = 
axis, the 0(3) FP, and the cubic FP for v > which should be the stable one. The stable 
cubic FP turns out to be very close to the 0(3) one, and the critical exponents are not 
distinguishable with the best 0(3) estimates. Indeed, FT estimates of their differences 
give [M] Uc — ^0(3) = —0.0003(3) and rjc — ''7o(3) = —0.0001(1). The rj conjecture would 
require rjc — ?7o{3) > 0; the above FT estimate of the difference r/c — ?7o(3)! although 
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favouring a negative sign, is not sufficiently precise to conclude that the t] conjecture 
fails. In conclusion, also in this case the rj conjecture is substantially consistent with 
the available results: for = 3 ?7cubic ~ ^0(3) > ''7ising > 0. 

The case N > 3. For > 3 the analysis of the six- loop series reported in Ref. [31] 
is consistent with the conjecture. One finds r/cubic ~ "/Rising, but the precision is not 
sufficient to determine which one is larger. 

The case N —>■ 00. For —>■ 00, keeping Nu and v fixed, one can derive exact 
expressions for the exponents at the cubic FP. Indeed, for ^ 00 the system can be 
reinterpreted as a constrained Ising model, leading to a Fisher renormalization [35] of 
the Ising critical exponents. One has 

V = Vi + 0{l/N), jj=^^ + 0{l/N), (19) 

where r//, uj, and aj are the critical exponents of the Ising model. Again the rj conjecture 
does not fail. 

3.3. Spin- density-wave model 

The spin-density-wave model is a rather complicate $^ model with five quartic 
parameters, which could describe the SDW-SC-to-SC phase transition in high-T^ 
superconductors (cuprates) [36l|. Its Hamiltonian density is 

n = \^,<^^\' + \d,<!>2\' + KI-KiP + |<f2p) + ^(1$ il' + \'^2\') 

where $a are complex A^-component vectors. The RG flow of this model has been 
investigated in Ref. [36]. The physical interesting cases are those for A^ = 2,3. The 
analysis reported in Ref. [5B] suggests the existence of a stable FP in both cases, with 
rather large values of 77: 77 = 0.12(1) and r] = 0.18(3) respectively (from MZM). There 
are also several unstable FPs in the RG flow, but all of them have smaller values of 
Tj. Therefore, the available d = 3 results of this spin- density- wave model support the rj 
conjecture. 

3.4. U{N) (g) U{N)- symmetric $^ models 
Let us now consider the Hamiltonian density 

'H = ii (9^$^) (9^$) + rtr <l>"f$ + ^ (tr ^ + ^tr (^t^) ^ (20) 

where $ij is a A^ x A^ complex matrix. The symmetry is U{N)l ® f/(A^)_R. In the 
case is also symmetric the sjTiimetry is U{N). This model has been introduced and 
studied [371 EH] because it is relevant for the finite-temperature transition in QCD, that 
is, the theory of the strong interactions. Indeed, for vq > the ground state leads to 
the symmetry-breaking pattern U{N)l U{N)r — > U{N)v (corresponding to QCD if 
the U{1)a anomaly is neglected). 
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For = 1, the model reduces to the 0(2)-symmetric theory. For > 3 no 
stable FP is found [38] . 

We focus on the case N = 2. One can easily identify two FP's. One is the Gaussian 
FP for u = V = 0, which is always unstable. Since for vq = the Hamiltonian becomes 
equivalent to the one of the 0(8)-symmetric model, the corresponding 0(8) FP must 
exist on the v = axis for u > 0. The 0(8) FP is also unstable because the f-term 
in the Hamiltonian represents a spin-4 perturbation with respect to the 0(8) FP, and 
such perturbations are relevant for any 0{K) models with K > 3 [15j. Within one-loop 
e-expansion calculations, no other FP is found [37J. By contrast, three-dimensional 
analysis of both MZM and MS expansions show the presence of a stable FP in three 
dimensions [38j- The corresponding value of rj is rather large: rj ^ 0.1, significantly 
larger than r]o{8) ^ 0.03. Therefore, it supports the r] conjecture. 

3.5. SU (A) -symmetric model 

Let us now consider the Hamiltonian density 

n = tT (9^$^) (9^$) + rtr <l>^<l> + ^ (tr ^ 

+ ^tr ^ + Wo (det <l>^ + det $) , (21) 

where $ is a complex and symmetric 4x4 matrix field. The symmetry of this model 
is S'f/(4). For vq > —^\wo\, the theory describes the symmetry-breaking pattern 
SU{4:) — * 5*0(4) [38], which is the appropriate symmetry-breaking pattern to describe 
transitions in a QCD-like theory with quarks in the adjoint representation. 

Close to four dimensions there are only two FPs, the Gaussian and the 0(20) FPs, 
which are both unstable. They remain unstable even at lower dimensions and thus are 
of no relevance for the critical behaviour. The three-dimensional RG flow has been 
investigated by field-theoretical methods based on perturbative approaches, within the 
MZM and MS schemes [38] • They show the presence of a stable three-dimensional 
fixed point characterized by the symmetry-breaking pattern SU(4) SO (4). The 
corresponding value of t] is rather larger: t] ^ 0.2, much larger than the values of t] of 
the Gaussian and 0(20) FPs which is 770(20) ~ 0.013. 

4. $^ theories of multicritical behaviours 

In this section, we discuss the extension of the rj conjecture to <l>^ theories describing 
multicritical behaviours, characterized by more than one independent correlation 
lengths. In this situation, the e-expansion indicates that the stable FP should be the 
one with the largest value of the trace of the rj matrix. 
We discuss this point within the $^ theory 

n = {djif + {0^2^ + ri0i' + r202' + )' + «2(02 )' + 02 (22) 

where 0i^2 are two 0(ni) and 0(^2) order parameters, with rzi and N2 real components 
respectively. The symmetry is 0(ni) ©0(^2). This theory describes the multicritical 
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behaviour arising from the competition of orderings with symmetries 0(ni) and 0(n2), 
at the point where the corresponding transition hnes meet in the phase diagram |39j . 
Such multicritical points arise in several physical contexts, for instance in anisotropic 
antiferromagnets, in high-Tc superconductors, etc.... See, for example, Ref. \15\ and 
references therein. 

The multicritical behaviour is determined by the RG flow in the quartic-coupling 
space when ri 2 are tuned to their critical values. Four FPs are found: the Gaussian FP, 
the isotropic 0{ni + FP (describing an effective enlargement of the symmetry), a 
decoupled 0{ni)-0{n2) FP (which describes effectively decoupled order parameters), 
and a biconal FP. The main properties of the three-dimensional RG flow are the 
following. 

(i) For ni + n2 > 4 the decoupled FP is stable. This can be inferred from non- 
perturbative arguments [ID] that show that the RG dimension of the perturbation 
Pw = 0102 th^t couples the two order parameters is 

= - + - - 3, (23) 

Z/i Z/2 

where i^i,z/2 are the correlation-length exponents of the 0(ni) and 0(^2) models. 
Inserting the numbers reported in Table [1] one finds y^ < for ui + n2 > 4 and 
any Ui, n2 > 1. 

(ii) Field-theoretical methods based on perturbative expansions (six-loop in the 
MZM and 0{e^) in the ^-expansion) show that the isotropic 0{ni + ^2) FP is unstable 
for ni + ^2 > 3 [15]. Therefore, only in the case of two Ising order parameters can the 
symmetry be effectively enlarged from Z2 © Z2 to 0(2), at the multicritical point where 
the Ising lines meet. 

(iii) For rii = 1 (Ising), n2 = 2 {XY), 0{6^) calculations show that the stable FP 
is the biconal FP [15]. Its critical exponents turn out to be very close to the 0(3) ones, 
in fact they are not distinguishable within the errors of the best estimates of the 0(3) 
critical exponents, see Table [H 

If we compare the isotropic and the decoupled FPs, the conjecture on the trace of 
the rj matrix should give 

niT]0(ni) + ^2^0(n2) > ('^l + ^2)^0(ni+n2) (24) 

for rii + n2 > 4, and 

VXY > msing (25) 
Moreover, from the point (iii) above, we should have 

tr r/biconal > ^?Ising + '^VXY (26) 

All these relations are verified by, or when the precision is not sufficient are consistent 
with, the best estimates of the exponent rj for the 0{N) models, see Tabled! 
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5. The Gross— Neveu— Yukawa model 

In this final section we discuss the Gross-Neveu- Yukawa (GNY) model ^1], and show 
that the infrared stable FP of its RG flow is the one characterized by the fastest decay 
of the critical two-point function of the boson field. 
The Lagrangian of the GNY model is 

^ = -J2 Ml A + 9o(T)i^r + ^[(^M^)' + m'a'] + ^uoa\ (27) 

i 

where ipi are Nf fermionic fields, and a is a real scalar field. The relation between the 
GNY and the standard Gross-Neveu model is discussed in Ref. |41j . 

Its RG flow can be investigated within the e-expansion. The RG functions at 
one-loop order are |5T] 

— f- 



i3u = -eu + — { -u" + Nug^ - 6Ng^ ] , (28) 



2 , ^ + 6 4 



where N = Nftil (the trace is in the 7-matrix space) is the total number of fermion 
components. In four dimensions trJ = 4, and thus = ANf in Eq. fl28|) . u and g are 
the MS renormalized couplings associated with uq and go, respectively. Note that at 
one-loop the RG (3 function for a completely general Yukawa coupling [12] derives also 
for a potential (see also [3]). 

The j3 functions of the GNY model have three FPs. Beside the unstable Gaussian 
FP, there is an Ising FP at 

= -^e, gl = 0, (29) 

which is also unstable. The infrared stable FP of the theory is the Gross-Neveu FP at 

^ 384iV7r2 2 ^ 167r2 

(Ar + 6)[(Ar-6) + (Ar2 + i32Ar + 36)V2]^' ^* N + Q^' ^ ^ 

We can now compare the corresponding RG dimensions of the scalar field, and check if 
the infrared stable FP is the one with the largest value of rj^j. Calculations of the scalar 
field RG functions show that the critical exponent 77^ is maximum at the GN FP: 

while ri„ = 0{e'^) at the Ising FP. 

The GNY model is soluble in the large- A^j limit for any d [41j, for a review see 
Ref. I30l. In this limit one finds 



V. = A-d + Oil/Nf), (32) 

and therefore ?7o- = 1 + 0{1/Nf) in three dimensions. The values of rjcr at the Ising FP 
are definitely smaller, for example r^^ ~ 0.036 in three dimensions, see Table [H and 
r/o- = 1/4 in two dimensions. 
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In conclusions, the above analytical results suggest that in the GNY model the 
infrared stable FP is the one that corresponds to the fastest decay of correlations of the 
scalar field, like in the $^ theories. 

Appendix A. Some proofs within the ^-expansion 

In the framework of the e-expansion, we prove two consequences of the property of 
gradient flow discussed in Sec. [2j (i) there exists at most one stable fixed point; (ii) the 
stable fixed point corresponds to the lowest value of the potential. Indeed, let us assume 
the existence of two fixed points corresponding to the parameters g* and g'*. We then 
consider the parameters g of the form 

g{s) = sg* + {l-s)g'\ < s < 1, (A.l) 

and the corresponding potential u{s) = U{g{s)). As the explicit form shows, at 
leading order u{s) is a third degree polynomial in s. The derivative 

u\s) = T.9'ais)^ = E(^?: - C)^ = E(^: - cmgis)) (a.2) 

vanishes due to the fixed point conditions at s = and s = 1: 'u'(O) = u'{l) = 0. Since 
u'{s) is a second degree polynomial, it then has necessarily the form 

u'{s) = As{l-s). (A.3) 

The second derivative u"{s) is given in terms of the matrix of second partial derivatives 
of U and, thus, the partial derivatives of the /9-functions, by 

u"{s) = E(^: - a'al^^^^ig; - g'C) = A{i - 2s) (A.4) 

In particular, for s = and s = 1 

At a stable fixed point, the matrix U" of partial second derivatives of U is positive. 
Thus, if g* and g'* are stable fixed points, A and —A are both given by the expectation 
value of a positive matrix and thus are both positive, which is contradictory: the two 
fixed points cannot both be stable. 

More generally, the sign of A characterizes, in some sense, the relative stability of 
these two fixed points. Let us assume, for example, A < which is consistent with 
the assumption that g* is stable. Then u'{s) < in [0, 1] and U{g{s)) is a decreasing 
function. Thus, 

f/(/) < U{g'*). (A.7) 

In particular, if g* is a stable fixed point, it corresponds, among all fixed points, to the 
lowest value of the potential. 
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